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ABSTRACT 

Shocks of supernova reniiiaiits (SNRs) are important (and perhaps the dominant) agents for produc- 
tion of the Galactic cosmic rays. Recent 7-ray observations of several SNRs have made this case more 
compelling. However, these broadband high-energy measurements also reveal a variety of spectral 
shape demanding more comprehensive modeling of emissions from SNRs. According to the locally 
observed fluxes of cosmic ray protons and electrons, the electron-to-proton number ratio is known to 
be about 1%. Assuming such a ratio is universal for all SNRs and identical spectral shape for all 
kinds of accelerated particles, we propose a unified model that ascribes the distinct 7-ray spectra of 
different SNRs to variations of the medium density and the spectral difference between cosmic ray 
electrons and protons observed at Earth to transport effects. For low density environments, the 7-ray 
emission is inverse-Compton dominated. For high density environments like systems of high-energy 
particles interacting with molecular clouds, the 7-ray emission is 7r°-decay dominated. The model 
predicts a hadronic origin of 7-ray emission from very old remnants interacting mostly with molecular 
clouds and a leptonic origin for intermediate age remnants whose shocks propagate in a low density 
environment created by their progenitors via e.g., strong stellar winds. These results can be regarded 
as evidence in support of the SNR-origin of the Galactic cosmic rays. 

Subject headings: radiation mechanism: non-thermal — gamma rays: ISM — ISM: supernova rem- 
nants — cosmic rays 



L INTRODUCTION 

Nearly a century after the discovery of cosmic rays 
(CRs), their origin remains one of the biggest funda- 
mental questions in the fields of high energy physics 
and astrophysics. It is widely believed that the super- 
nova remnants (SNRs) are one of the most probable 
candidate sour ces of the Galactic CRs below the so - 
called "knee" (Bhat ct al.' ToSS"; 'Erlvkin & WolfcndaL 
l2001al|bl E [2002a; Hillas 2005; Katz & Waxman 200 
However, there are currently no observations that can 
directly verify such a conjecture. 

The multi-wavelength observations of SNRs, espe- 
cially the high energy 7-rays from the ground-based 
atmospheric Cerenkov telescope arrays and space-based 
telescopes, provide powerful tools to probe the particle 
acceleration mechanism in SNRs (lEnom oto et al.l 120021: 
Aharonian et al.ll200l: lUchivama ct al. 20071 : iTianet all 
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20071 120081 I2OI0I: iNeronov fc Semikoz,,2012E . However, 
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whether the nature of these 7-ray emission from SNRs 
is predominantly hadronic or leptonic is still a ma tter of 
debate (e. g., [A haronian ct al. 2 006a; Ber ezhko fc ^ Volkl 



20M iButt et al. 200 8: iLiu et al. 2008 ; iPlagal 
Morlino et 31112009; .Fang et al.l .2009; .Fan et al.l 
Ellison et al.l 120101 : Zirakashvili fc AharonianI 
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Berezhko fc Volkl l2010l : lYuan et al.l l2011al) . That is to 



say, SNRs are known particle accelerators, but it is not 
clear whether they dominate the observed Galactic CR 
flux, especially for nuclei, on the Earth. 

The observational 7-ray spectra of SNRs also show 
significant diversity. Recent Fermi observations of the 
young SNR RX J1713. 7-3946 shows a very hard spec- 
trum in GeV energy range, 1.50 ±0.11, which implies an 



inver se-C ompton (IC) origin of the 7-rays (jAbdo et al.l 
I2OIII see llnoue et al.l ()2012t ) for an alternative explana- 
tion). On the other hand, for all of the SNRs inter- 
acting with molecular clouds (MCs), the GeV-TeV 7- 
ray spectra are generally very soft and seem to better 
agree with the 7r°-decay model, although the model of 
bre msstrahlung emission from electron s can not be ruled 
out JAbdo et al.ll2010dl. ^009l [2010lbl[an : lAiello et al.l 
I2OI2I : iGiuliani et al.ll2010l ). It is natural to ask whether 
there is a common understanding of these 7-ray signa- 
tures of the SNRs. 

In this paper we illustrate that the locally observed 
CRs and the 7-ray emission of SNRs can be naturally 
understood in a unified picture. Assuming both the 
protonqj and electrons are accelerated in SNRs and the 
SNRs are the dominant sources of Galactic CRs, one can 
derive the electron-to-proton ratio at the source based 
on the locally observed spectra of protons and electrons. 
Then the 7-ray emission of SNRs will depend primar- 
ily on the environmental gas density, and the observed 
diversity of 7-ray spectra of SNRs can be attributed to 
variations of this density. 

2. ELECTRON-PROTON RATIO AT THE SOURCE 

The observational flux of protons at ^GeV is about two 
orders of magnitude higher than that of electrons, which 
implies the electron-proton ratio at the source should 
be of the order 1% since the energy loss of electrons is 
negligible in this low energy range. This result is well 
known (e.g., iCohen fc Ramatvl 119731 : lApparao fc Daniell 

^ The heavier nuclei which may play a similar but less important 
role than protons for 7-ray emission are not discussed here. 



[19771 ILevinsOT][l99l lErivkin fc WolfCT^dali [2002bh . In 
this section we derive the injection parameters of the 
Galactic CRs at the sources by reproducing the observed 
spectra on the Earth, considering the detailed propaga- 
tion model. 

After production at sources, charged energetic parti- 
cles propagate diffusively in the random magnetic field 
of the Galaxy. The overall convection and reaccelera- 
tion due to scattering with the random magnetohydrody- 
namic waves may also change the distribution function of 
CR particles. Furthermore, the interactions between the 
CRs and the gas, interstellar radiation field, and mag- 
netic field, will lead to fragmentation, catastrophic or 
continuous energy losses of these particles. The trans- 
port of C Rs from the sou rces to the Earth is generally 
complex (|Butt et al.ll2008l see lStronget all ([2007) for a 
recent review of the CR propagation) . 

In this work we limit our study to CR protons 
and electrons only. We a,dopt the GALPROP codjl 
(jStrong fc Moskalenkoll998l:lMoskalenko fc StrondlQQSl) 
version v50p to calculate the propagation of the CRs. 
The diffusion-reacceleration frame, without convection, 
is assumed. The main propagation parameters are Do — 
6.59 X 10^*^ cm^s-i, S = 0.30, va = 39.2 km s~^ and 
Zh = 3.9 kpc, which are derived through t he fit to the 
B/C i°Be/^Be, Carbon and Oxygen data (jTrotta et al.l 
1201 ID . The spatial distribution of the CR sources in the 
cylindrical coordinates is 



exp I --^ 



f{R,z)(x l—j exp — 

(1) 
where axis symmetry has been assumed and Rq = 8.5 
kpc is the distance of solar system from the Galactic 
center, Zs « 0.2 kpc is the scale height of the source 
distribution, a = 1.25 and b = 3.56. Such a source func- 
tion is similar to the SNR spatial distribution but tuned 
base d on the Fermi obs ervations of diffusive Galactic 7- 
rays (jTrotta et al.ll2011[ ). 

The injection spectral shape as a function of momen- 
tum p (or rigidity) is assumed to be a broken power-law 
function 
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which is identical for both electrons and protonfQ. The 
ratio of the normalization between electrons and protons, 
usually called K^p, is taken as a free parameter. The bro- 
ken power-law inj ection spectrum is required to fit the 
observed CR da ta (jStrong et al.ll2004l:lTrotta et al.ll2011l: 
Liu et al.]l2012[) as wel l as the 7-ray da t a (lAbdo et al 



2010c 



l2010ellal H lAiello et al.ll2012t INeronov et ah 



2012f l. It was proposed that strong ion- neutral collisions 
near the shock front may lead to the spe ctral break of 
accelerated particles around ^ 10 GeV (jMalkov et al.l 
120051 12011[ ). Alternatively the escape effect of particles 
from/into finite-size region may also give a break at sev- 
eral GeV (jOhira et al.l(20TTI : iLi fc Chenll201(t [Ml . 

Compared with the PAMELA observations of the 
CR proton (jAdriani et al.l l2011bf ) and electron spectra 

^ http://galprop.stanford.edu/ 

^ For the sake of simpHcity, we ignore the spectral evolution 
and potential spectral difference b etween high energy electrons and 
protons escaping from the SNRs IjErlvkin fc Wolfendale|[2002bl '). 




Fig. 1. — The expected fluxes of CR protons and electrons at 
the Earth, for the same spectral shape of the inje cted particles, 
compa red with the PAMELA observational data llAdriani ct al.| 
l2011ai rblV We adopt two parameter settings to calculate the elec- 
tron spectrum: for solid line the magnetic field is the canonical 
one adopted in GALPROP and Kep ~ 1.3%; for dashed line the 



magnetic field is two times larger and K„ 



1.9%. 



(jAdriani et al.l l2011al) . we find ai = 1.80, 02 = 2.52 
and PbrC = 6 GeV can give an acceptable fit to the 
data, where c is the speed of lighl[3, as shown in Fig. [TJ 
For energies below ~ 30 GeV, the sol ar modulation with 
force -field approximation is employed (jGleeson fc AxfordI 
|1968[ ). To fit both the proton and electron spectra with 
the same injection spectrum, we need different modu- 
lation potentials, as labeled in Fig. [T] This may be 
due to the rest-mass or sign of charge dependence of 
the modulation effect (|Clem et al.lll996[) . To reproduce 
the absolute fluxes, we find the electron-to-proton ra- 
tio Kep is about 1.3%. Such a result could be rcpro- 
duc ed with the numerica l simula tion (jLevinsonI [1994) . 
See lErlvkin fc Wolfendald (|2002bf ) for more possible ex- 
planations. 

A better fit to the electron data can be obtained 
through increasing the magnetic fielc^ by a factor of 2 
(correspondingly K^p « 1.9%), as shown by the red- 
dashed line. Even with an identical spectral shape 
at injection, the electron spectrum is much softer 
than that of protons due to their energy loss in the 
transport processes from the source regions to Earth, 
which is quite differe n t from the scenario explored by 
lErlvkin fc Wolfendali (l2002bf ). where SNRs produce an 
electron distribution softer than the proton distribution. 

As a comparison, the independent fit to the proton 
and electron data gives ai = l-91i ctg = 2.40, phrC = 10 
GeV for protons (jTrotta et al.ll20lH) . and a i = 1.50 
02== 2.56, PhrC = 3.60 GeV for electrons (jLiu et al.l 
I2OI2I ). Given uncertainties of the propagation parame- 
ters, astrophysical inputs, solar modulation etc, we con- 
sider these spectral fits to be consistent with each other. 

* We are not dedicated to discuss the spectral hardening of CR 
nuclei above ~ 200 GeV reported by ATIC/CREAM/PA MELA 
llPanov et al.l [20071 : \Khn et al.|[20TOl : lAdriani et al.|[20n bl). which 
may im ply the superposition of different sour ce spectra ( Yuan ct ahl 
2011d) or a nearby new component of CRs I IErlvkin fc Wolfeudalg 



2013 : fl'houdam fc Horandcl 20121V 



^ Note the chang e of magnetic field may affect the synchrotron 
radiation JStrong et al.ll2011l ) . What we employ here is an example 
to include the uncertainties of the propagation model. The full dis- 
cussion of a self- consistent propagation model is beyond the scope 
the present work. 



3. GAMMA-RAY EMISSION OF SUPERNOVA 
REMNANTS 

We now investigate the 7-ray emission of SNRs adopt- 
ing Kep ^1% and the spectral parameters at the source 
as derived in the previous section. In general there 
are three major components of the 7-ray emission: IC 
and bremsstrahlung radiation by electrons and the 7r°- 
decay emission by protons. Therefore we further need 
the knowledge about the radiation background and the 
environmental gas density. For the sake of simplicity, we 
assume IC scattering with the cosmic microwave back- 
ground radiation only in this work. The scattering with 
infrared and o ptical lig ht may make the IC 7-ray spec- 
trum broader (|Porter et al. 2006), but is not expected to 
affect the qualitative discussion here, at least for those 
far away from the Galactic center. Then the only param- 
eter determining the SNR 7-ray spectra will be the gas 
density. 

We classify the SNRs into three groups with low, 
medium, and high ambient gas densities. The lack of 
thermal X-ray emission of RX J1713. 7-3946 gives an 
upper limit of the g as den sity of 0.02 — 0.03 cm^'^ 
(|Cassam-Chenai' et al.l l2004t lYuan et all l2011al) . For 
RX J0852.0-4622 ASCA X-ray data implies a gas den- 
sity n < 0.03((i/lkpc)-i/2/-i/2 cm'^, with d the dis- 
tance and / the fi lling factor of X-ray emitting volume 
(|Slane et al. 1120011) . Although the derivation of the upper 
limit relies on assumptions such as the ionization equilib- 
rium and gas temperature, we take them as the typical 
examples of low density SNRs. Adopting n = 0.01 cm~^, 
together with the CR spectra given by Equation (2) and 
Kep ~ 1%, we show the calculated 7-ray spectra of these 
two SNRs in Fig. [21 The overall normalization of the 
7-ray luminosity is determined by the observational flux 
of each source. To be consistent with the high-energy 
spectral cutoff behavior of both SNRs, we employ an ex- 
ponential cutoff term with Ec ~ 60 TeV of the electron 
spectra. The cutoff might be due to the balance of accel- 
eration and the cooling in the vicinity of the SNR. For 
protons the cutoff could be higher and is assumed not 
to enter the energy range discussed here. A remarkable 
signature of the 7-ray spectrum is that, the spectrum is 
very hard and the luminosity of TeV emission is higher 
than that of the GeV emission. The 7-ray emission of 
these SNRs is IC dominated. 

Some relatively young SNRs have a moderate gas den- 
sity, such as Cassiopeia A and Tycho. The average gas 
densit y of Cassiopeia A i s estimated to be about 4.4 
cm~^ () Arava fc Cuil [20T0I) . For Tycho an upper limit 
n < 0.6 cm~'^ was derived from the absence of thermal 
X-ray emission from the brig ht outer rim of the remnant 
(jCassam-Chenai et al.|[2007D . The density in the inner 
region of the remnant can be much higher. Here we take 
these two SNRs as examples of medium density sample 
and adopt n = 1 cm~'^ 0. The calculated spectra are 
shown in Fig. [3] We can see that for this kind of sources 
the GeV emission is 7r°-decay dominated and the TeV 
emission is IC dominated. The luminosities in GeV and 
TeV bands are comparable in this case. 

® Although these sources c learly show com plex source structure 
with multiple emission zones I IAtovan &: Derr ncr 2012), in this pa- 
per we still consider the simple one zone emission model to capture 
the dominant features. 



Finally we discuss the case with a high density, typi- 
cally for the SNR-MC interacting systems. The gas den- 
sity in the MCs can easily reach 10^ — 10^ cm~^. As 
an illustration we adopt n = 100 cm^'' in this study. 
The expected 7-ray spectra together with the GeV- TeV 
observational data of eight SNR-MC systems are shown 
in Fig. m In this case the GeV- TeV 7-ray emission is 
TT^-decay dominated, and the 7-ray spectrum is generally 
very soft. 

We may need to estimate the impact of energy losses 
in the dense clouds on the CR proton spectrum. For n — 
100 cm~'^ ISM, the ionization en ergy loss time scale of 
proto ns is about 10^ yr at 1 GeV (jStrong fc Moskalenkcl 
Il998f ). and the pion-production energy loss time scale is 
about 10^ yr. For most of these SNRs the ages are es- 
timated to be less than several tens kyr, so we expect 
that for the energy range interested in this work (GeV- 
TeV) the energy losses of the protons are not important. 
For the low energy particles (several tens MeV) and high 
energy electrons (> TeV), the energy losses may be im- 
portant and need to be considered in future modeling. 
As for the contribution to 7-rays from the secondary e , 
it is only important for E^ < 10 MeV compared with the 
j pjp induced Tr'^-decay component, even for very old SNR 
(IFang fc ZhandlMl . 

In Fig. [5] we show the relation between the photon in- 
dex r and the gas density n, for the SNRs studied in this 
work. The parameters of the SNRs are compiled in Table 
[TJ The photon index F is fitted using the observational 
data between 1 GeV and 1 TeV, with single power-law 
function. For the SNR-MC interacting system whose gas 
density is not well known we assume a value of 10^ — 10^ 
cm""^. A trend showing the correlation between F and 
n can be seen from Fig. \5\ We also show the theoreti- 
cal expected result based on the unified model with the 
solid line. The model is consistent with the observational 
data. 

It is encouraging that the results in Figs. [2M] show 
rough consistence with the observational data, in support 
of our relatively simple interpretation of the complicated 
7-ray spectral behaviors of SNRs. Note here we have 
not tried to precisely fit the observational data because 
variations of the spectral and env ironmental parameters 
are expected for different sources (jFerrand fc MarcowithI 
120101: lYuan et al . 2011b3). With slight adjustment of 
these sou rce parameters w e can easily get better fit to 
the data (|Li fc Chenll2"oTl) . 

4. CONCLUSION AND DISCUSSION 

In this work we propose a unified model to explain the 
7-ray emission of SNRs. In the model, by assuming that 
SNRs produce identical high energy electron and pro- 
ton spectral shape and are the major sources of the low 
energy CRs (below the "knee"), the electron-to-proton 
number ratio Kep at the sources is derived to be ~ 1% 
according to a realistic CR propagation model described 
with the GALPROP code. With such a K^p value we 
calculate the expected 7-ray spectra for various SNRs 
with different environmental parameters. Qualitatively 
the observed diversity of 7-ray spectra of different SNRs 
can be naturally understood with different gas densities. 
For low density environments the 7-ray emission is IC 
dominated, while for high density environments the 7- 
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Fig. 2.— Expected 7-ray spectra for SNRs RX J1713.7-3946 (left) and RX J0852. 0-4622 (right). The gas density is adopted to be 
n = 0.01 cm-3. Refe rences of the observ ational d ata — RX J1713.7-3946 : Fermi HAbdo et al.|[20Til) . HESS IjAharonian et al.ll2007bh : RX 
J0852.0-4622: Fermi UTanaka et al.ll20Tll '). HESS HAharonian et al.ll2007g ). 
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Fig. 3. — Same as Fig. [2] but for Cassiopeia A (left) and Tycho (right). The gas density is a dopted to be n = 1 cm •^. Refe rences of 
the ob servational data — Ca ssiopeia A: Fermi llAbdo et al.ll2010bl 'l. MAGIC l| Albert et al.ll2007bl 1. VERITAS HAcciari et al.|[20TOl ): Tycho: 
Fermi l IGiordano et al.ll20Tl ). VERITAS l lAcciari et al.ll2011l~ 



TABLE 1 

Parameters of the SNRs: Name, R.A., Dec, distange, age, gas density, photon index (between 1 GeV and 1 TeV) and 

REFERENCES. 



Name 



R.A. 



Dec. 



d (kpc) Age (kyr) n (cm '') 



Ref. 



RX J1713.7-3946 


17'^ 13"" 50= 


-39°45' 


1.0 


1.6 


<0.03 


1.77 ±0.08 


1,2 


RX J0852.0-4622 


08''52™00= 


-46° 20' 


0.75 


1.7-4.3 


<0.03 


1.89 ±0.04 


3,4 


Cassiopeia A 


23''23™26= 


-l-58°48' 


3.4 


0.32 


4.4 


2.19 ±0.04 


5-8 


Tycho 


00''25™18= 


-1-64° 09' 


2.5-3.0 


0.44 


<0.6 


2.18 ±0.11 


9-12 


W28 


18''00™30= 


-23°26' 


2.0 


35 - 150 


— 


2.67 ±0.03 


13,14 


W41 


18''34™45= 


-08°48' 


3.9-4.5 


60 - 200 


— 


2.27 ±0.03 


15,16 


W49B 


ighiirUQ^s 


-l-09°06' 


8-12 


1-4 


— 


2.73 ±0.05 


17,18 


W51C 


19''23™50'' 


-|-14°06' 


6.0 


30 


— 


2.58 ±0.04 


19-21 


IC443 


06'' 17^00= 


-1-22° 34' 


0.7-2.0 


3-30 


60 - 240 


2.61 ±0.04 


22-24 


CTB 37A 


17''14™06= 


-38°32' 


6.3-9.5 


— 


— 


2.48 ±0.07 


25-27 


G8. 7-0.1 


18''05™30= 


-21°26' 


4.8-6.0 


15-28 


— 


2.23 ±0.05 


28,29 


G359.1-0.5 


nh^^m^QS 


-29° 57' 


7.6 


>10 


— 


2.63 ±0.02 


30,31 




Abdo et al.| fl2011|): 




et al.| 201 


l7b|):(3)|Tanakaetal.| 


20111): (4)|Aharonianet al 


1 2007a| 
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ray emission is Tr^-dccay dominated. The model predicts 
that 7-ray spectra in low density environments are gen- 
eral harder than those of SNR-MC interaction systems. 



Since strong thermal emission is expected from shocked 
dense media, we expect relatively weaker thermal emis- 
sion from remnants with harder 7-ray spectra than those 
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Fig. 4. — Same as Fig. [2] but for S NR-MC interacting systems. The gas density is ad opted to be n = 100 cm ^. References of 
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Fig. 5. — The photon index F (between 1 GcV and 1 TeV) versus 
the gas density n of the 12 SNRs studied in this work. The solid 
Une is the model expected result. 

with softer 7-ray spectra. Such a simple, self-consistent 
model, if further validated by observations, supports the 
SNR-origiu of the low-energy CRs. 

The age could also be a parameter affecting the 7-ray 
emission of SNRs, and it can be coupled with the density 
parameter. For example, considering the progenitors, the 
density is in the intermediate range for very young rem- 
nants, lowest f or middle-age SN Rs, and highest for very 
old remnants (jDwarkada J 120051 ) . An improved unified 
model of course needs to consider more factors, includ- 
ing the effect of multiple emission zon es as observed in 
some sources (jAtovan fc Dermeill2012f) . which will affect 
the detailed spectral fit. 

In general, for Kf,p '^1%, the hadronic component will 
always dominate the bremsstrahlung component, which 
is distinguishable from the model th at bremsstrahlung 
may dominate the 7-ray emission (jAtovan fc Dermerl 
l2012f l. and could be tested by the observation of 7-rays 
in lower energy range (< 100 MeV). 

When calculating the 7-ray emission of the SNRs, we 
employ the particle spectra which are the same as the 
injection spectra giving rise to the locally observed ones. 
However, we should keep in mind that the CR spectra 
accelerated by the source at specific epoch may be differ- 
ent from that injected into the Mil ky Way, which shoul d 
be the time integrated spectrum (|Caprioli et al.|[2010D . 
A more detailed modeling may need to tak e into account 
the evolution hi story of the SNR s (e.g . . iSturner et al. 
19971 : [Rerooldi [20081: iLeeet al.l mm iFang fc Zhang 



200a iFinke fc Dermeij 120121 : lYan et d\\\2Ql'£ . Further 
more, the electron-to-proton ratio is assumed to be a con 
stant in this work, independent of the sources. Such an 
assumption can break the degeneracy between medium 
density n and K^p ■ It is a strong assumption but is not in 
conflict with observations. Detailed study of individual 
sources is needed to verify this assumption. 

The magnetic field of the shocked emission region is 
an important parameter but docs not appear explicitly 
in our model. However, for the multi-wavelength mod- 
eling including the synchrotron emission component, the 
average magnetic field might be well-determined via a 



detailed spectral fit ()Fan et al.ll2010D . Indeed with the 
increase of the magnetic field, the IC component will 
be suppressed for a given synchrotron fiux due to the 
decrease of the number of energetic electrons. The gas 
density (or number of CR protons and therefore 1/Kep) 
needs to be increased in this case to account for the ob- 
served 7-ray flux. The 7-ray emission is likely domi- 
nated by pionic emission due to higher gas density or 
proton fiux in this case, and it will be difficult to de- 
tect the IC component. Another effect for the high mag- 
netic field is the cooling of accelerated electrons via syn- 
chrotron emission. Thus the electron spectrum should be 
in contradiction with our assumption of a unified spec- 
tral shape. Except for this particular case our study will 
remain valid as far as this strong magnetic field does not 
dominate the overall particle acceleration in SNRs. In 
cases where IC dominates the 7-ray emission, the mag- 
netic field can be directly derived through the radio to 
X-ray emission of the SNRs. The magnetic field is usually 
weak (~ lOfxG) with an energy density only a factor of 
a few times higher than that of the background photons, 
according t o the sy nchrotro n X-ray to IC 7-ray lumi- 
nosity ratio (jLiu et a l. 2008). iLagage fc Cesarskvl (|1983l ) 
showed that with such a weak field, SNRs can barely ac- 
celerate charged particles up to 10^^ cV. Although the 
discussion in this work may still hold {E < 100 TeV) 
even in the case of weak magnetic field, for CRs to be 
accelerated up to the knee via the diffusive shock accel- 
eration, magnetic fields need to be amplified furth er by 
some mechanisms (e.g.. lBel]|[200l iGuo et al.ll2012| ). In- 
homogeneity of t he magnetic field and mu lti-zone accel- 
eration scenario (jAtovan fc Dermeil 120121 ) may alleviate 
challenges to the weak field case. Better determination of 
the magnetic field, together with the gas density param- 
eter, will provide crucial tests of the assumptions made 
in this work. 

Finally, the radio observations of SNRs indicate that 
the electro n spectrum is about E~'^ with a remarkable 
dispersion ('Cl ark fc Cas"^ill [19761 IBoedan et all [HH; 
[Schlickciser 20Q2), which is not exactly the same as that 
inferred in Sec. 2. We also expect dispersion of other 
parameters characterizing the energetic particle distri- 
bution, which may be used to improve the fit to the ob- 
served 7-ray spectra. The present work gives a zero-order 
approximation to the problem of the SNR 7-ray emission 
and the origin of CRs. Some common features of the 7- 
ray emission of different population of SNRs are revealed. 
Further works about the details may be helpful to refine 
the present model. 
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